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ABSTRACT 

The recently discovered transitional millisecond pulsar system J1023+0038 exposes a crucial evo- 
lutionary phase of recycled neutron stars for multiwavelength study. The system, comprising the 
neutron star itself, its stellar companion, and the surrounding medium, is visible across the electro- 
magnetic spectrum from the radio to X-ray /gamma-ray regimes and offers insight into the recycling 
phase of millisecond pulsar evolution. Here, we report on multiple-epoch astrometric observations 
with the Very Long Baseline Array (VLBA) which give a system parallax of 0.731 ± 0.022 milliarc- 
seconds (mas) and a proper motion of 17.98 ± 0.05 mas yr _1 . By combining our results with previous 
optical observations, we are able to use the parallax distance of 1368lgg pc to estimate the mass of 
the pulsar as 1.71 ± 0.16 Mq, and we are also able to measure the 3D space velocity of the system as 
126 ± 5 km s _1 . Despite the precise nature of the VLBA measurements, the remaining ^3% distance 
uncertainty dominates the 0.16 Mq error on our mass estimate. 

Subject headings: Astrometry — pulsars: individual(J1023+0038) — techniques: interferometric — 
pulsars: general 



1. INTRODUCTION 

"Recycled" millisecond pulsars attain their high spin 
frequ encies by accreting matter from a companion donor 
star (lAlpar et al.l ll982: Bhattac harva fc van den Heuvell 
1991). The mass-transfer phase is identified by variable 
X-ray emission from the heated gas in the pulsar accre- 
tion disk, and this phase is thought to be relatively well 
characterized theoretically. However, the end of accre- 
tion and birth of the radio-visible millisecond pulsar is 
still poorly understood. The discovery of a millisecond 
pulsar system whose optical variability suggests it had 
an acc retion disk in 2001 (PSR J1023+0q3l& hereafter 
J1023: lArehibald et al.l 12001 IWang et al.l I2009T) offers a 
unique opportunity to study this short-lived transitional 
regime. 

The J1023 system comprises a 1.69-ms radio pulsar 
with a nondegenerate companion in a 4.75-hour orbit. 
It displays unusual observational characteristics which 
are consistent with its apparent evolutionary stage - the 
pulsar dispersion measure varies on timescales of min- 
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utes to weeks, and long-term variation in the orbital 
period is also apparent. Eclipses of the pulsar emis- 
sion are observed, usu ally at particular orbital phases 
(| Archibald et al.ll2009l ). The influence of material lib- 
erated from the companion object and interacting with 
the pulsar wind is likely responsible for all of these phe- 
nomena. The J1023 system, then, offers a rich source of 
information for the study of binary evolution and inter- 
acting wind physics. 

As a compact radio source of ~mJy brightness, J1023 
is a suitable target for astrometric observations using 
very long baseline interferometry (VLBI), which can 
provide extremely high distance accuracy via the mea- 
surement of annual geometric parallax. Pulsar astro- 
metric campaigns with the Very Long Baseline Array 
(VLB A) in the US ([Chatterjee et al.l [20091: [Brisken et all 
2002) and the Lon g Baseline Array (LBA) in Australia 
([Deller et~a l. 2009) have obtained parallax measurements 
with accuracies at the level of 20 /zas. 

2. OBSERVATIONS AND DATA PROCESSING 

2.1. VLBI observations 

Due to the very high resolution of VLBI observations 
(<10 mas), even a detection of moderate significance 
yields a very precise position centroid, and systematic 
errors relating to calibr ation typically domina te the as- 
trometric error budget ([Chatterjee et al.ll2004n . By far, 
the dominant contribution to systematic calibration er- 
rors at frequencies below several GHz is the unmodeled 
differential ionosphere between the calibrator and target 
sources. This contribution can be reduced by observ- 
ing at higher frequencies or utilizing a calibrator very 
close to the target on the sky. As Table Q] illustrates, 
J 1023 has an average (although highly variable) flux of 
several mJy at 1.6 GHz. Our 1.6 GHz astrometric obser- 
vations detected J1023 with a significance ranging from 
7 to 63, and since J1023 possesses a typical steep pulsar 
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TABLE 1 
Observation details 



Observation 


Observation 


MJD 


Pulsar flux 


group 




dates 






(mjy) a 


2008a 


Dec. 


16 


54816 




2.0 


2009a 


Feb. 


14 


54876 




0.7 


2009b 


May 


17, Jun. 12 


54968, 


54995 


3.0, 1.1 


2009c 


Nov. 


21, Dec. 7 


55156, 


55172 


0.9, 1.9 


2010a 


May 


15, 29 


54816, 


55346 


0.4, 3.1 


2010b 


Nov. 


17, Dec. 3 


55517, 


55533 


2.5, 1.1 



a Averaged over the pulsar period 



spectrum, higher frequency observations would not have 
regularly detected the pulsar. Accordingly, minimizing 
the ionospheric error was only possible by minimizing 
the calibrator-target separation. 

To this end, we made use of an "in-beam" calibrator 
(whose angular separation from the target is less than 
the size of the telescope primary beam, and can thus 
be observed simultaneously with the target). By ana- 
lyzing archival Very Large Array observations we iden- 
tified the flat spectrum object PMN J1023+0024 as a 
likely compact source, and confirmed its suitability with 
exploratory VLBA observations. PMN J1023+0024 is 
bright (peak flux density 80 mJy/beam at 1.6 GHz) 
and almost unresolved (deconvolved size 2.5x1.5 mas). 
It is separated from J1023 by just 15' (J2000 position 
10:23:38.77064 +00:24:16.6690; estimated uncertainty 
0.5 mas in each coordinate). 

We conducted 10 astrometric observations of J1023 
spread over a 2-year period from late 2008 to late 2010. 
In each observation, both right and left circular polar- 
ization were recorded in four observing bands of width 8 
MHz, resulting in a total data rate of 256 Mbps/antenna. 
The bands were placed adjacent to one another and 
spanned the frequency range 1651.49 - 1683.49 MHz. 
After two exploratory observations, the final eight obser- 
vations were clustered in four pairs, with each pair being 
observed at a time close to parallax extrema to improve 
our sensitivity to the parallax signature. The observa- 
tions are summarized in Table [1] The observations were 
phase referenced to the bright VLBA calibrator J1024- 
0052, and the pointing center for the target source was 
placed midway between J1023 and PMN J1023+0024. 
Since the in-beam calibrator was relatively bright, an 
unusually long cycle time (10 minutes) back to the pri- 
mary calibrator was used, as this "phase referencing" 
was used solely to improve the absolute position of the 
in-beam calibrator and to refine the amplitude calibra- 
tion. A total of 180 minutes on-source time was obtained 
per epoch - after the excision of data during the pulsar 
eclipse (discussed below) , the resultant continuum image 
sensitivity was typically 75 /^Jy, although a number of 
epochs suffered instrument failure at one or more anten- 
nas, slightly reducing sensitivity. 

2.2. Supporting timing observations 

In contrast to the very stable timing behavior exhib- 
ited by most millisecond pulsars, J1023 exhibits small, 
apparently random fluctuations in its orbital period on 
short (week to month) timesca les, making extrapola - 
tion of timing solutions difficult ((Archibald ct al. 2009). 
In order to predict the pulsar rotational phase accu- 



rately and employ pulsar gating to boost S/N (see e.g., 
IDeller et al.l [20091. we carried out regular timing obser- 
vations of J1023 in support of our VLBA program. These 
observations are briefly summarized below. 

For the firs t VLBA epochs we used the timing solution 
published in [Archibald et al.l ([2009P . For later epochs, 
we observed J1023 for an hour every few weeks with the 
Arecibo t elescope and the ASP coherent dedispersion 
backend ([Demoresti [20071: iFerdmanl [2001 at 1.4 GHz. 
This yielded adequate ephemeridcs until early 2010, 
when damage to the Arecibo telescope was discovered, 
requiring zenith angle restrictions that precluded obser- 
vations of J1023. J1023's erratic dispersion measure vari- 
ations made phase connection to ongoing 350 MHz Green 
Bank Telescope observations unreliable, and so for the 
final epochs we made use of data from the J1023 tim- 
ing programs carried out using t he Lovell Telescope at 
Jodrell Bank (jHobbs et al.l 120041 ) and the Westerbork 
Synth esis Radio Telescope (WSRT; iKaruppusamv et all 
l2008lh 

An additional complication is that J1023 is un- 
detectable at 1.4 GH z for about 30% of the orbit 
([Archibald et a.1.1 12009T ) . After correlation, we discarded 
all integrations that fell between orbital phases 0.15 and 
0.45 (around 60 minutes of on-target time per observa- 
tion) using a simple or bital ephemeris that assumes a 
constant orbital period ([Archibald et al.ll20Toh . 

2.3. Correlation and data processing 

From the first three observation groups, all epochs ex- 
cept 2009 May 17 were correlated using the VLBA hard- 
ware correlator. The epoch at 2009 May 17 and the final 
6 epoch s were correlated using the DiFX software cor- 
relator ([Deller et al.l l201lh . In each case, two correlator 
passes were employed; the first without pulsar gating at 
the position of PMN J1023+0024, and the second using 
the pulsar gating scheme described below at the position 
of J1023. An averaging time of 2 seconds was used for 
both passes. 

Pulsar gating involves zeroing the correlator output 
when the pulse phase is outside a specified window, pre- 
venting the addition of unnecessary noise into the visibil- 
ity outputs. J1023 has a wide, double peaked pulse shape 
([Archibald et al.ll2009f ). which limits the gain from pulse 
gating to a factor of ~1.7 when a simple on/off gate is 
used. The erratic timing behavior of J1023 precludes the 
use o f more advanced "matched filtering" (IDeller et al.1 
120071 ). which could otherwise have boosted the gate gain 
to a factor of ~2.2. 

The visibility data produced by the correlator were 
reduced using AIPS^I. utilizi ng standard scripts based 
on the ParselTongue package ([Kettenis et al.l [20061 ). Af- 
ter loading the data and flagging known bad data (in- 
cluding time ranges affected by the eclipse), the vis- 
ibility amplitudes were calibrated, fringe fitting and 
bandpass calibration were performed using the bright 
source J0927+3902, and the delay and amplitude so- 
lutions were further refined with fringe fitting and am- 
plitude self-calibration using the out-of-beam reference 
source J1024-0052. Subsequently, phase-only correc- 
tions were generated from the in-beam calibrator PMN 
J1023+0024. This and all future steps were undertaken 

10 http://www.aips.nrao.edu/ 
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on the Stokes I data. For all reference sources, a com- 
bined model was formed based on the data from all 
epochs; these models were not permitted to vary between 
epochs. 

Finally, the J1023 data were split and averaged in fre- 
quency, leaving a single visibility point per band for each 
integration, and the visibil ity data were w ritten to disk 
and imaged using difmap ()Shepherdl 119971 ) with natural 
weighting (which maximizes the image sensitivity). The 
resultant images (one per 8-MHz band, Stokes I, here- 
after the "single-band" images) were fitted with a single 
gaussian component using the ALPS task JMFIT, as was 
a combined image formed using all subbands (hereafter 
the "combined" image). For two epochs (MJD 54876 
and 55331) the pulsar was not detected in the single- 
band images, and only the combined image was formed 
and fitted for position. 

2.4. Astrometric data reduction 

The astrometric error in the position measurement for 
each epoch is principally composed of a random compo- 
nent (the error in fitting a position centroid in the image 
plane) and a systematic component (from the unmodeled 
delay errors at each telescope due to the atmosphere and 
ionosphere). The pulsar brightness varies by an order of 
magnitude throughout our 10 observing epochs; random 
errors dominate the error budget when it is faint, whilst 
for the remaining epochs the systematic error dominates. 
This fact, along with the highly variable atmospheric and 
ionospheric conditions between epochs, means that the 
distribution of astrometric errors in right ascension and 
declination will not be Gaussian. 

A straightforward fit to the 10 positions obtained from 
the combined images yields a reduced \ 2 value of 2.2, 
indicating that the errors on a whole are underesti- 
mated due to the unaccounted-for systematic contribu- 
tion. Simply inflating the error budget by scaling the po- 
sition errors at each epoch to obtain a reduced \ 2 value 
of 1.0 is inappropriate, since that assumes the system- 
atic error contribution is correlated with the random er- 
ror contribution at each epoch. A better, although still 
imperfect approach is to add an equal error in quadra- 
ture to all epochs. In order to obtain a reduced \ 2 value 
of 1.0, the required error magnitude was 50 fi&s, added 
solely to the right ascension coordinatcEE This astromet- 
ric error corresponds to an average residual path length 
error of 2 mm. At our observing frequency of 1.6 GHz, 
the ionospheric correction will be by far the dominant 
source of path length error, as current ionospheric mod- 
els sample the sky coarsely in time and area. During our 
observations, the differential ionospheric delay correction 
between the calibrator and the target ranged from to 
80 picoseconds (0 to 25 mm), with typical values of ^20 
picoseconds (6 mm). A path length error of 2 mm, then, 
corresponds to an average error of around 30% in the 
ionospheric correction. The fitted parameters obtained 
when a constant systematic error contribution of 50 fias 
is added to each epoch (fit "A") are shown in the first 
data column of Table [5] The position errors shown in 
parentheses for right ascension and declination are for- 
mal errors shown for comparative purposes only. The ac- 

1 The fit varies insignificantly if a larger total error of 70 fj,as is 
divided equally between right ascension and declination. 



TABLE 2 

Fitted astrometric parameters for J1023. 



Parameter 



Fit A (combined) 51 Fit B (bootstrap) 



Right ascension (J2000) b 


10:23:47.687198(2) 


10:23:47.687196(2) 


Declination (J2000) b 


00:38:40.84551(4) 


00:38:40.84551(4) 


Position epoch (MJD) 


55000 


55000 


/i a (mas yr _1 ) 


4.76 ± 0.03 


4.79 ± 0.04 


W (mas yr^ 1 ) 


-17.34 ± 0.04 


-17.35 ± 0.05 


Parallax (mas) 


0.731 ± 0.022 


0.730 ± 0.028 


Distance (pc) 


1368 ttl 


1370 t\\ 


i>T (km s _1 ) 


11711 


U7±l 



a Values from the combined fit (fit A) are used in the analysis. 
b The errors quoted here are from the astrometric fit only and 
do not include the ~0.5 mas position uncertainty transferred 
from the inbeam calibrator's absolute position. 

tual errors in the absolute position coordinates are dom- 
inated by the absolute position of the reference source 
PMN J1023+0038 and are estimated at 0.5 mas. 

An alternate method to estimating the astromet- 
ric errors is to use bootstrapping to estimate the ac- 
tual form of the error distribution without imposing 
any prior constraints. B ootstrapping is described in 
lEfron fc Tibshiranl (|1991|) and has been previously used 
in the estim ation of errors in pulsa r astrometry obsery - 
ables (e.g., IChatteriee et all 120091 : IDeller et al.l 1200911 . 
In essence, it involves creating a large number of test 
datasets, where each dataset is constructed by sampling 
with replacement from the pool of measured astrometric 
positions. The astrometric observables are fitted once 
from each test dataset and a large sample of tests is used 
to build a histogram of the fitted values for each observ- 
able. The results of a bootstrap error analysis for J1023 
using positions from the single-band images (fit "B" ) are 
shown in the second data column of Table [2J 

The disadvantage of the bootstrapping technique is 
that a relatively large number of observed positions are 
required to avoid small-sample-size uncertainty. This 
necessitates using the single-band position fits for the 
epochs where they are available, but the relatively low 
significance of the detections becomes problematic. The 
positions (and associated errors) obtained from image 
plane fits become increasingly unreliable at low S/N ra- 
tios; in fact, the weighted mean position obtained by 
averaging the single-band positions can differ from the 
combined image position fit by up to 0.5tr. Accordingly, 
we would expect the bootstrap error values to be more 
conservative than those in fit A, and this is indeed seen, 
although the difference is small and the fitted values are 
consistent. Hereafter, we use the values and errors from 
fit A in all further analysis. The motion of J1023 over 2 
years and the parallax signature after the subtraction of 
the best-fit proper motion are plotted in Figure [T] 

3. IMPLICATIONS OF THE ASTROMETRIC RESULTS 

The measured parallax of 0.731 ± 0.022 /ia.s yields a 
distance to the J1023 system of 1368lgg pc. The very 
high precision of the measur ement means that th e Lutz- 
Kelker bias in the distance (TVerbiest et al.|[20Tfj|) is neg- 
ligible for this system. Combined with the total mea- 
sured proper motion of 17.98 ± 0.05 mas yr^ 1 , the trans- 
verse velocity is calculated as 117 ± 4 km s _1 . The 
distance is over twice that predicted by the NE2001 elec- 
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Best astro metric fit 

O Error ellipses for observations 



2 4 6 

Right Ascension offset (mas) 



,-rI 



Best astrometric fit 
Measured positions 



800 54900 55000 55100 55200 55300 55400 55500 55600 
Time (MJD) 



Best astrometric fit 
Measured positions 



800 54900 55000 55100 55200 55300 55400 55500 55600 
Time (MJD) 



Fig. 1.— The astrometric fits to the J1023 VLBI data. 
Top panel: motion relative to the reference position of 
10:23:47.687198+00:38:40.84551 over the two years of observations, 
with the best fit shown by the dashed line. Middle panel: The par- 
allax signature in right ascension (proper motion subtracted). Bot- 
tom panel: The parallax signature in declination (proper motion 
subtracted). 



tron density distribution model (600 pc; iCordes fc Laziol 
2002). Accordingly, the average electron density along 
the line of sight is lower than expected, just 0.010 
cm" 3 . The systemic radial velocit y for J1023 reported by 
iThorstensen fc Armstrong! (|2005f ) is negligible (1 ± 3 km 
s _1 ), and so the 3D system velocity is equal to the trans- 



verse velocity calculated here. Correction for peculiar 
solar motion and Galactic rotation using a flat rotation 
curve and the current IAU recommended rotation con- 
stants (Ro = 8.5 kpc, 0n = 220kms _1 ) alters this value 
slightly to 126 ± 5 km s _1 . 

Having a precise distance and velocity makes it pos- 
sible to calculate corrections to the observed pulsar 
spin period derivativ e P a pp due to the Shklovskii ef- 
fect (IShklovskiil 119701) and accelerat ion in the Galactic 
potential (e.g., iNice fe Tavlorl 11995 s ) . After these cor- 
rections are subtracted to obtain the true spin-period 
derivative P, this can be used to obtain accurate val- 
ues for quantities such as the spin-down luminosity E, 
the magnetic field at the neutron star surface B sur f, and 
the characteristic age r c . Using the supporting timing 
observations made with the Lovell Telescope, we obtain 
P ap p = (6.83 ±0.05) x 10~ 21 . The contribution due to 
the Shklovskii effect is Pshk/P = M 2 D/c, where /i is the 
proper motion, D is the pulsar distance and c is the speed 
of light. Substituting the distance and velocity derived 
above, and taking P as 1.6879874 ms (fArchibal d et al.l 

l2009h . we obtain P shk = (1.81 ± 0.04) x 10" 21 . Finally, 
we calculate the net effect of acceleration in the Galac- 
tic gravitational potential on P ap p to b e (—2.8 ± 0.3) x 
10~ 22 , following INice fe Tavlorl (119951) and taking the 
vertic al acceleration component from lrlolmberg fe Flvnnl 
(|200l . We used R = 8.5 kpc and 9 = 220 km s" 1 , and 
assumed 10% errors in these constants and the K z force 
law of iHolmberg fe Flvnnl (|2004l ). Thus, the intrinsic P 
is (5.30 ±0.07) x 10" 21 . 

This revised P is consi derably lower than the prelim- 
inary value presented in lArchibald et al.1 J2009), with 
consequent effects for E, B sur f, and r c . E is reduced to 
4.3 x 10 34 ergs s" 1 , B surf is reduced to 9.6 x 10 7 G, and 
t c is increased to 5.0 x 10 9 years. The value for B sur f is 
noteworthy given that less than 10% of recycled pulsars 
have B surf < 10 s G. 

3.1. The mass of J 1023 

When IThorstensen fe Armstrong! (|2005f ) carried out 
optical modelling of the J1023 system, they noted that 
the mass of the companion could be inferred from the 
distance, assuming that it fills its Roche lobe. Essen- 
tially, the distance combined with the multiband pho- 
tometry determine the temperature and size of the com- 
panion, and the assumption links this to the size of the 
Roche lobe and therefore the mass of the com panion M c . 
Specifically, IThorstensen fc Armstrong! ((2005) fit models 
including Roche geometry and companion temperature 
distributions to the observed multiband light curves and 
find that the distance to the system (under the assump- 
tion of Roche-lobe-filling) is: 



1/3 

D = 2.20 [ I kpc 



M, 



Table 5 of IThorstensen fc Armstrong! (|2005D lists 49 
sets of model parameters that fit the observed optical 
data. Testing the above relationship on all 49 reveals 
a scatter of only 1% on the coefficient (2.20 kpc) over 
the broad range of parameters they consider plausible 
enough to fit. 
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The pulsar timing measurements of lArchibald et al.l 
( 2009) showed that the system mass ratio Mpsr/M c = 
7.1 ± 0.1. Combining these two relationships, we can ex- 
press the mass of the neutron star in terms of the distance 
to the system as: 



Mpsr - (7.1 ±0.1) x 



D 



2.20 ± 0.02 kpc 



M, 



©• 



Taking our distance measurement of 1368^3g pc and 
adding all three uncertainties in quadrature, we can 
obtain Mpsr = 1.71 ± 0.16M Q , where the distance 
uncertainty remains the dominant contributor to the 
error budget. This result is consistent with recent 
analyses which have showed that recy cled pulsars have 
masses of approxi mately 1.5 ± O.2M (jOzel et al.ll2012t 
IZhang et alJl201ll) . 

The assumption of Roche lobe filling by the companion 
of J1023 remains a concern. While it seems likely that 
the active episode ex perienced by the J102 3 system in 
2001 (summarized in lArchibald et akl I2009D was due to 
Roche lobe overflow, we have regrettably limited infor- 
mation about the nature of that episode. It is possible 
that the active episode was actually the result of some 
other form of mass transfer, perhaps due to very strong 
winds from the companion, though we deem this unlikely. 
If the 2001 activity did indeed involve Roche lobe over- 
flow, iThorstensen fc Armstrong] (|2005D point out that 
the Kelvin-Helmholtz relaxation time for the compan- 
ion is much longer than the several years since the active 
episode, implying that the companion should still fill its 
Roche lobe. Regardless, if the companion's radius is a 
fraction f (necessarily less than one, or excess material 
would fall away) of the Roche lobe radius, the pulsar 



mass is approximately (1.71 ± O.16M0)f -3 . That is, if 
the companion substantially underfills its Roche lobe, the 
mass estimate we present is a lower limit for the pulsar 
mass. 

It would be valuable to have observations that 
tested whether the comp a nion fills its Roche lobe. 
IThorstensen fc Armstrong] (|2005l ) found that models 
with a heated Roche-lobe-filling companion fit the opti- 
cal light curves very well, but they did not explicitly test 
models that did not fill the Roche lobe. Sensitive optical 
photometric observations might reveal more structure in 
the light curves, but phase-resolved high-quality spec- 
troscopy is an even more promising avenue for probing 
the companion's size and shape. 
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